Novel Fermi liquid behaviour of atomic mixtures 
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We show that the measurement of the quasiparticle lifetime can reveal new Fermi liquid regimes in atom 
mixtures with density and/or mass imbalance. In these recently created mixtures, the existence of two Fermi 
energies leads to novel and unusual Fermi liquids where the quasiparticle decay rate is no longer quadratic in 
energy and temperature. For certain densities and mass ratios, the decay rate exhibits a linear scaling similar to 
that of marginal Femi liquids, whereas in other regimes it exhibits a plateau. It is demonstrated that this novel 
Fermi liquid behaviour can extend to rather high temperatures which are within current experimental reach. The 
novel Fermi liquid regimes can be detected with radio-frequency spectroscopy. 



Introduction. — Landau's concept of the Fermi liquid [1] 
has played a crucial role in our understanding of many sys- 
tems ranging from liquid helium, electrons in metals to atomic 
nuclei and quark-gluon plasmas. In ultracold atomic Fermi 
mixtures, several studies have already measured Fermi liquid 
thermodynamic [2, 3] and spectroscopic properties [4-7]. A 
hallmark of conventional Fermi liquid physics is the quadratic 
dependence of the quasiparticle decay rate on excitation en- 
ergy and temperature. Here we show that cold atom gases 
with density and/or mass imbalance can realise new Fermi liq- 
uid regimes, where the decay rate scales linearly or is even a 
constant as a function of energy and temperature. The linear 
scaling is reminiscent of the so-called marginal Fermi liquid 
- a phase which plays an important role in the investigation 
of unconventional and high T c superconductivity. For these 
reasons we argue that the measurement of the quasiparticle 
lifetime would be a new and important experimental result 
since, not only would it show the conventional quadratic be- 
haviour which has never been observed in atomic gases, but it 
would also reveal the presence of new unconventional regimes 
in these mixtures. 

The reason for the appearance of new scaling laws is that 
in a system with two types of atom (which we will denote by 
"t" and "|") of different masses and/or densities, the Fermi 
energies of the two species are in general unequal leading 
to the existence of two energy (and temperature) scales. We 
will assume here that £p± <C £p<\. The difference can be a 
result of a large density imbalance n^ <C n^ with equal mass 
or large mass imbalance ^> with equal density. This 
leads to the existence of a new intermediate regime in energy 
£ fi %t wm ch increases in size with mass or density 

imbalance. While the quadratic behaviour of conventional 
Fermi liquids is due to the effect of Pauli blocking on both 
species participating the scattering, the key property of this 
new intermediate regime is that the Pauli blocking of the 4- 
atoms is unimportant so that the phase space for interspecies 
scattering depends only on the t atoms and this leads to novel 
Fermi liquid behaviour. We will show that, for the 
equal mass mixture, the Pauli blocking of t atoms simply re- 
sults in a linear decay rate in temperature while for am|»/n|, 
equal density mixture, energy conservation restricts the phase 
space for scattering of the t atoms to a narrow region around 



the Fermi surface, making Pauli blocking of the t" atoms irrel- 
evant, leading to a temperature-independent decay rate. Ex- 
perimentally, Fermi mixtures with density imbalance [2, 8] or 
mass imbalance [9-11] have been created which raises the 
prospect of observing the effects we discuss here in the near 
future using for instance radio-frequency (RF) spectroscopy. 

Model. — We consider a homogeneous gas of two species of 
fermions denoted a =t,| with masses < and densities 
— n l> f rom which we define the Fermi momenta kp G = 
(6n 2 n G ) 1 ^ . The key quantity we study in this paper is the 
decay rate 1 jx p of the t quasiparticles with momentum p and 
energy £ p ^. To lowest order in the scattering processes, the 
decay rate can be written as 

= 2nU 2 £ £ 8 (£ pt + £ K - £ p _ qt - £ k+q; ) 

Z P k q 

X[ftk;(l -^k+q|)(l -^p-qt) + ( 1 -^k^k+q^p-qt] C 1 ) 

where n kG = (e^ ko + 1) _1 is the Fermi function. We have 
defined £ pcy = £ pcy — jl G with \i G the chemical potential, and 
P = 1/7 with T the temperature (we set ks = h = 1). The 
energy is taken to be £ pG = p 2 /2m G where m G is the effective 
mass. The parameter U is the effective interaction between 
the t and the | atoms, and we have neglected interactions be- 
tween identical atoms for simplicity. In the strong coupling 
regime, one can extract the value of the effective mass and 
U from Monte Carlo, variational and thermodynamic argu- 
ments [12-14]. As is usual in Fermi liquid theory, we will 
assume that U is only weakly dependent on energy and mo- 
mentum below £p^. In the weak coupling regime, we have 
U = 27la/m r where m r = m^m^/(m^ + m^) is the reduced 
mass and a the scattering length for the interaction between 
the two atom species [15]. 

For analytic investigation, it is convenient to rewrite (1) in 
terms of the imaginary part of the Lindhard function of the 4- 
atoms which is defined as [14] 

f d 3 k 

Im^(^,co) = -7r J j^(n H -n k + qi )S(cQ-e k + qi + £ H ). 

(2) 

To do this, we first recast the 5 function in (1) in the 
form <5(£ pt + £ki ~ £ P -qt - ^k+ q ;) = J^dco8(co - e pt + 



2 



£ p _ q ^)5(o) — £k+q; + £14). We also use the Fermi func- 
tion identities ^(1 -ft k+q ;) = (n K -n k+q ;)/(l -e~^) and 

(1 -%;H+q; = -^k+q|)/(^ w - 1), with (0 = e k+q; - 
£14. Finally, the angular integral over q is f£l q S(co — £ p | + 
£ p-qt) = 2nmt/pq with —pq/m^ — q 2 /2m^ < CO < pq/m^ — 
q 2 /2m^ and we obtain 



m t \U\ 2 
2n 2 p 



POO p(0 + 

/ dqq dCDlmXi(q,CD)F(CD,£tf,^) 
JO J co- 



(3) 

where F(to,£ pt ,^ t ) = (1 + £? P( fl >-5 P t))-i(i - e^)" 1 + (1 + 
e -P(co-^)yl^ e pa) _ and = ± pq / mr - ( f/2m r . 

Before proceeding, let us briefly discuss the difference be- 
tween the situation considered here and that of conventional 
Fermi liquids. In a conventional Fermi liquid, the low energy 
condition % p = £ p — ji <C £f and T <C £^ ensures that one can 
use lmXi(q, co) co/q in (3) leading to quadratic scaling of 
the decay rate with energy and temperature. Here, the Lind- 
hard function of the | atoms has a different behaviour in the 
intermediate regime, which will result in different power laws. 

Zero temperature. — In the following we will assume that 
t; P t > without loss of generality since T p (£ ) is an even func- 
tion for |£ I <C £pt' m m * s case ' me back-scattering term [the 
second term in (3)] vanishes at zero temperature. The region 
of integration in (3) is determined by three conditions: from 
the Bose factors we have < CO < % p y, CO < C0+(q) and < 
q <2p. Using the zero-temperature expression for the Lind- 
hard function [16] lmxi(q, co) = -m 2 £ F ^/4nq [0(1 -v 2 _)(l - 
v 2 ) -0(1 -v 2 )(1 -v 2 )] with v±=m i co/qk n ±q/2k n , we 



obtain 



1 \U\ 2 m^m 2 £ F i r^p rtrf, co + 



87T 3 /? 



/ dq 

Jo Jo 



dco 



x[0(l-v 2 _)(l-v 2 _)-0(l-v 2 + )(l-v 2 + )]. 



(4) 



(i) The low energy regime <C £ F ± <C £f\- When the 
energy is small compared to both Fermi energies, we can use 
Im;^(g, co) = —m 2 co/47lq which gives 



\U\ 



'm^mfkfi 

87T 3 p 



(5) 



This is the usual quadratic dependence of the decay rate on ex- 
citation energy characteristic of a conventional Fermi liquid. 
Indeed, when = we recover the well-known expression 
for the damping rate of a quasiparticle at T = [16]. 

(ii) The intermediate regime £p^ < ^ < J£ F ^ <C £f\ with 
y=4(k Ft /k n - l)(^ t /^ n + m t /m;)/(l+m t /m;) 2 where 
J£ F ^ is defined as the CO coordinate of the intersection of the 
&>+ and v_ = — 1 curves (note that for this region to exist we 
must have 7 > 1 which sets a condition on the density and 
mass ratios). We obtain 



\U\ 2 m^m\n Fi 
2%p 



(6) 



The linear dependence of the decay rate on excitation en- 
ergy is also characteristic of marginal Fermi liquids [17]. In 
Fig. (1) we plot the decay at zero temperature as a func- 
tion of the excitation energy. We have chosen the param- 
eters m^/m^ — 113/6 corresponding to a mixture of 173 Yb 
and 6 Li atoms [10], and kp^/kpi = 2. These parameters give 
£ Ft/ £ Fi = 115 corresponding to a large regime of interme- 
diate energies. Both the usual Fermi liquid (quadratic) and 
marginal scalings (linear) are clearly visible in Fig. (1). 
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Figure 1: (color online). Zero temperature decay rate \/t p (units of 
\U\ 2 m^m 2 £ p ^/87t 3 ) of an t quasiparticles as a function of excitation 
energy with m^/m^ = 173/6 and kpt/kp^ = 2. The black solid 
curve is the result of a numerical integration of (3) and the red dashed 
curves are given by (5) and (6). 

Non-zero temperature, degenerate case. — In this regime, 
we have < T <C 7>| and we consider the case ^ = 0. From 
the thermal distribution functions in (3) it follows that the inte- 
grand is significant only in the range \ co\ <T. Since T 
we can approximate (3) by 



1 



Ft 



\U\ 2 mt 

27l 2 k 



Ft 



r 2k Ft r 

Jo qdq L 



2ImXi(q,C0)dCQ 
(ePa + lXl-e-P®)' 



(7) 



(i) T <C Tp^. In this case, we can use the zero temperature, 
low-energy expression Im;^(g, co) = —m 2 co/47lq to obtain 



\U\ 2 mtm 2 k Fi ^ 2 



T kpt %7lkpt 



(8) 



Here the system has conventional Fermi liquid behaviour and 
we recover the standard result when Tp± = Tpt [16]. 

(ii) Tp\, «K< Tpt- The gas of | atoms is now clas- 
sical, and we can use the expression Im;^(g, co, T) = 

- nn± y/2m$ j% e-™ 2 ™^/ 2 * 2 e^P^l sinh(/3 CO /2) /q [14]. 
The decay rate can then be written as 



1 



1 



with 



'(<) 



dx- 



-x 2 /4y 2 

coshx 



(9) 



(10) 
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We now discuss two important special cases of (9)-(10). First 
we consider the case of a highly polarized system of two spin 
states of the same atom, i.e. = = m and <C n^. In 



t 2 ) for t <C 1 we find that 



this case, we have ^Jk F ^/2mT = yj 7>| / f ^> 1. Using /(oo) 
V^ln2, we obtain 



\U\ z m z n\ 
:21n2^ ±T. 



(11) 



■Ft 



The linear temperature behaviour is again reminiscent of a 
marginal Fermi liquid. In Fig. 2, we plot the decay rate 1 / T^ Ft 
as a function of T for several density imbalances. In Fig. 2 
(a) we see the Fermi liquid low temperature behaviour. Note 
that the curves for kp^/kp^ = 3, 6, 10 overlap in this range. In 
Fig. 2 (b) we see the appearance of linear scaling as in the 
marginal Fermi liquid. As density imbalance increases, the 
range of temperatures exhibiting linear scaling increases. 

The second case we consider is that of an equal density 
mixture of heavy and light atoms (m^ <C and n = = n±) 



so that Jkl/lm^T = y/T Fi /T <C 1. Since I(t) ~ ^fnt 2 {\ 
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Figure 2: (color online) decay rate 1/T^ n in units of 
U 2 rrr > £^^kp^/S7i^kp^ as a function of temperature with — — 
m for different density imbalances kp^/kp^ = 1,3,6, 10. The black 
(dark) solid curves are the numerical integration of (3) while the red 
(light) solid curve in (a) and (b) is given by (8) and (11) respectively 



1 _ |£/| 2 m t ?2|& n 7>; 



(1 



(12) 
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which shows that the decay rate is constant to leading order in 
Tp^/T. This peculiar behaviour can be understood as follows. 
Since T <C Tp^, we have v t h^ <C Vp^ with v t h^ = y/TT Jrn\, 
and vp^ = kp^/m^, and the 4- atoms are moving very slowly 
compared with the t atoms. The decay rate is then domi- 
nated by the motion of the t quasiparticle: T _1 = n^G sc vp^ 
where where <7 SC is the scattering cross section. Using <7 SC = 
\U\ 2 m 2 /n for m; > [14], we get T _1 ~ \U\ 2 m^k F ^/7i, 
which is precisely the leading constant term of the above ex- 
pression. In Fig. 3, we plot the decay rate of an t excitation 
with momentum kp^ as a function of temperature. The emer- 
gence with increasing mass ratio m^/m^ of a plateau where the 
rate is independent of temperature is clearly visible in Fig. 3 
(a). Note that this plateau extends to temperatures well above 
7>| for large mass ratios even though the simple kinetic argu- 
ment breaks down above T F ^ since we then have to take into 
account the changes to the f Fermi surface which also affect 
the scattering rate. Instead, as we show below, the reason the 
plateau continues to higher temperatures is because of a can- 
cellation between the changes to the forward and backward 
scattering terms. 

Non-zero temperature, classical case (p = kp^). — In the 
regime Tp^ <C T, the t distribution is also classical. 

(i) 7>| «T« Tp^m^/m^ . Assuming that U continues 
to be weakly dependent on energy and momentum (as oc- 
curs for example when k F ^\a\ <C 1), the plateau will per- 
sist to a much larger temperature scale, given by v t h; ~ vp^ 
(i.e. T ~ Tp^m^/m^) as can be seen from the numerics 
(see Fig. (3)). This seems to indicate that the Pauli block- 
ing of the t atoms plays no role for the decay rate of an 
t quasiparticle since the plateau survives independently of 
whether the t atoms are degenerate (T <C 7>|) or classical 
(T ^> 7>|). To better understand this unusual behavior, we 
plot in Fig. 3 (b), the forward and backward scattering con- 
tributions to l/i£ Ft with m^/m^ = 100 [11]. We find that, 
when T <C Tp^, the forward and backward scatterings con- 
tribute equally to 1/t^, which can be understood by tak- 
ing £ pt = in (3). In the regime 7>t < T <C 7> t (m;/m t ), 
the backward scattering begins to decrease while the for- 
ward scattering increases, keeping however their sum con- 
stant. This shows that Pauli blocking indeed affects the for- 
ward and backward scatterings but not the sum of the two. 
This constant sum is due to the vanishing energy transfer in 
the decay process of the t quasiparticle with a large mass 
imbalance ( co max /T ~ max(m|/m^, y/Tp^/T) < 1). Thus 
by rewriting n H (l - n k+q ;)(l - n p _ qt ) + (1 - n H )n k + qi = 

(n H - % +q ;)/(^ - 1)^(1 - *p- q t) + flp-tf] * (n H - 
n k+qi) /(eP® — 1), we see that the t atoms play no role in the 
integrand of (3), explaining why the plateau survives indepen- 
dently of whether the t atoms are degenerate or classical, i.e., 
the Pauli blocking is irrelevant in this case. 
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(ii) T F ^m^/m^ <C T. Using the classical limit of the chemi- 
cal potential for fixed particle number in (3), i.e., [l/T —> —°° 
when T —> °o, we find 1 / T^ Ft is again given by (9), but with 

I(t)=r d ye-y 2 [^l^ (13) 

where t = ^Jk F ^/2m\T . When t <C 1 the integral can be eval- 
uated straightforwardly, which gives I(t) = 2tm 2 /m^m^ and 
the decay rate becomes 

1 2V2\U\ 2 m 2 r n\ r~ 
= ' ' r l Vf. (14) 

It is interesting to compare (14) with (12). First we note that 
the Fermi momentum k F ^ in (12) has been replaced by \fT 
in (14) as expected for the high temperature regime. Second, 
when the two expressions are equated we obtain that the cross- 
over between the intermediate T behaviour given by (12) and 
the high T behaviour given by (14) occurs for T ~ Tp^m^/m^ 
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Figure 3: (color online) decay rate in units of 

\U\ 2 k F m 2 /327t 3 (m^m^) 1 / 2 as a function of temperature for = 
obtained from a numerical integration of (3). The left figure shows 
the appearance of a plateau for m^^> m^. The right figure compares 
the numerical results (solid curves) with (8), (12), and (14) (dashed 
curves) in the low, intermediate and high temperature limits respec- 
tively. It also shows the contributions of the forward and backward 
scattering terms of (1). 



in agreement with the analysis above. Thus, the region where 
the damping rate is independent of temperature when = 
and ^> extends to temperatures much higher than 7>|. 
This effect is clearly illustrated in Fig. (3). Importantly, it 
makes the experimental observation of this new Fermi liquid 
regime significantly easier as it emerges already for fairly high 
temperatures, when the mass imbalance is large. 

Polaron case. — Given its experimental importance, we fi- 
nally calculate the collision rate of the minority particles |. 
When the density of n\ —> 0, i.e., in the polaron limit, we find 
from (1) 

^=T^' 2 ™.-K e ' 2+5 ?B (15) 

for £^ <C £^t an d T <C 7>|. This indicates that the po- 
laron shows normal Fermi liquid behavior, consistent with the 
literature[13]. 

Experimental probes. — We now discuss possible ways to 
detect the novel Fermi liquid regimes analysed above. Radio- 
frequency (RF) spectroscopy has proven to be a very success- 
ful method to probe the single particle properties in cold atom 
gases, and it should be well suited to detect the novel Fermi 
liquid regimes. First we note that the intermediate temperature 
regime Tp± «T« 7>| has already been achieved by the Inns- 
bruck group using a mixture of 6 Li and 40 K atoms [6]. This 
group furthermore probed the single particle properties of the 
I atoms ( 40 K) by coupling them to a non-interacting state us- 
ing RF spectroscopy. This was done both in the linear regime 
and in the highly non-linear regime. In the non-linear regime, 
the 40 K atoms performed Rabi oscillations between the inter- 
acting and the non-interacting state, and their single particle 
decay rate could be extracted from the damping of the oscil- 
lations. A similar experiment using RF spectroscopy on the 
t atoms should be able to detect their decay rate. In order to 
single out the new Fermi liquid regimes, one can block out the 
low lying states in the Fermi sea of t atoms from participating 
in the RF spectroscopy by coupling them to a non-interacting 
state with a filled Fermi sea with a Fermi energy slightly be- 
low the Fermi energy of the t atoms. The low lying states are 
then inert to the RF probe due to Fermi blocking of the non- 
interacting state. Alternatively, one could use momentum- 
resolved radio frequency spectroscopies which have also been 
implemented experimentally [6, 7, 18]. For solid state sys- 
tems, deviations from standard Fermi liquid behaviour of- 
ten show up in electrical resistivity, specific heat and mag- 
netic susceptibility measurements. For atomic gases, transport 
coefficients have been measured using spin transport experi- 
ments [19], and we expect these novel Fermi liquid properties 
to also show up in such experiments. 

Discussion and conclusions. — We have shown that the ex- 
istence of new energy and temperature scales in imbalanced 
systems can lead to novel Fermi liquid regimes. In particu- 
lar we have focus sed on the lifetime of the quasiparticles and 
shown that its measurement can reveal scaling laws which are 
not quadratic. In future work, we plan to address transport 
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properties such as the momentum relaxation which should 
also reflect the unusual characteristics of these regimes. 

Acknowledgements. — We would like to thank M. Baranov 
and P. Massignan for discussions. ZL and CL acknowledge 
support from the EPSRC through grant EP/I018514/1. 



* Electronic address: z.lan@soton.ac.uk 

[1] L. D. Landau, Sov. Phys. JETP 3, 920 (1957). 

[2] S. Nascimbene, N. Navon, K. J. Jiang, F. Chevy, and C. Sa- 
lomon, Nature, 465, 1057 (2010). 

[3] K. Van Houcke, F. Werner, E. Kozik, N. Prokofev, B. Svistunov, 
M. J. H. Ku, A. T. Sommer, L. W. Cheuk, A. Schirotzek and M. 
W. Zwierlein, Nat. Phys. 8, 366 (2012). 

[4] S. Nascimbene, N. Navon, S. Pilati, F. Chevy, S. Giorgini, A. 
Georges, and C. Salomon, Phys. Rev. Lett. 106, 215303 (2011). 

[5] A. Schirotzek, C.-H. Wu, A. Sommer, and M. W. Zwierlein, 
Phys. Rev. Lett. 102, 230402 (2009). 

[6] C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, P. Massig- 
nan, G. M. Bruun, F. Schreck and R. Grimm, Nature, 485, 615 
(2012). 

[7] M. Koschorreck, D. Pertot, E. Vogt, B. Frohlich, M. Feld and 

M. Kohl, Nature, 485, 619 (2012). 
[8] Y. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, Nature, 

451, 689 (2008). G. B. Partridge, W. Li, R. I. Kamar, Y. Liao, 

and R. G. Hulet, Science, 311, 503 (2006). 
[9] E. Wille, F. M. Spiegelhalder, G. Kerner, D. Naik, A. 

Trenkwalder, G. Hendl, F. Schreck, R. Grimm, T. G. Tiecke, 

J. T. M. Walraven, S. J. J. M. F. Kokkelmans, E. Tiesinga, and 



P. S. Julienne, Phys. Rev. Lett. 100, 053201 (2008). L. Costa, J. 
Brachmann, A.-C. Voigt, C. Hahn, M. Taglieber, T. W. Hansen, 
and K. Dieckmann, Phys. Rev. Lett. 105, 123201 (2010) 

[10] H. Hara, Y. Takasu, Y. Yamaoka, J. M. Doyle, and Y. Taka- 
hashi, Phys. Rev. Lett. 106, 205304 (2011). A. H. Hansen, A. 
Y. Khramov, W. H. Dowd, A. O. Jamison, B. P.-Swing, R. J. 
Roy, and S. Gupta, arXiv: 121 1.2267. 

[11] Even for equal mass spin mixtures, an unequal mass mixture 
can be created by e.g. imposing a spin dependent optical lattice 
which can change the effective mass of one of the species by 
extremely large values. 

[12] P. Massignan and G. M. Bruun, Eur. Phys. J. D 65, 83 (2011). 

[13] F. Chevy, Phys. Rev. A 74, 063628 (2006); C. Lobo, A. Re- 
cati, S. Giorgini, and S. Stringari, Phys. Rev. Lett. 97, 200403 
(2006); N. V. Prokofev and B. V. Svistunov, Phys. Rev. B 77, 
020408 (2008); G. M. Bruun and P. Massignan, Phys. Rev. Lett. 
105, 020403 (2010). 

[14] G. M. Bruun, A. Recati, C. J. Pethick, H. Smith, and S. 
Stringari, Phys. Rev. Lett. 100, 240406 (2008). 

[15] In the case of atomic gases we also assume that 3 -body recom- 
bination rates are small and do not affect the lifetime. 

[16] G. F. Giuliani and G. Vignale, Quantum theory of the electron 
liquid, Cambridge University Press, 2005. 

[17] CM. Varma, Z. Nussinov, W. van Saarloos, Phys. Rep. 361, 
267 (2002). Note however that, unlike in this reference, here 
in this paper we assume that perturbation theory is valid and 
Z=l + 0(a 2 ). 

[18] J. T. Stewart, J. P. Gaebler and D. S. Jin, Nature 454,744 (2008). 
[19] A. Sommer et al, Nature 472, 201 (2011); A. Sommer, M. Ku, 
and M. W. Zwierlein, New J. Phys. 13, 055009 (201 1). 



